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Quantum search algorithm tailored to clause satisfaction problems 
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Many important computer science problems can be reduced to clause satisfaction problem. We 
are given n Boolean variables Xk and m clauses Cj where each clause is a function of values of some 
of XkS. We want to find an assignment i of XkS for which all m clauses are satisfied. Let fj{i) be 
a binary function which is 1 if J**' clause is satisfied by the assignment i else fj (i) = 0. Then the 
solution is r for which f{i = r) = 1, where f{i) is the AND function of all fj{i)s. In quantum 
computing, Grovers algorithm can be used to find r. A crucial component of this algorithm is the 
selective phase inversion A of the solution state encoding r. A is implemented by computing f{i) 
for all i in superposition which requires computing AND of all m binary functions fj{i)s. Hence 
there must be coupling between the computation circuits for each fj{i)a. In this paper, we present 
an alternative quantum search algorithm which relaxes the requirement of such couplings. Hence it 
offers implementation advantages for clause satisfaction problems. 

PACS numbers: 03.67.Ac 


I. INTRODUCTION 

Grovers algorithm is used to search an item satisfying 
certain properties out of a database of N items [l|. Let 
the index i {i G {0, 1, ..., A^ — 1}) denote these items. 
Consider a quantum system of n = log 2 N qubits with 
Hilbert space of dimension N, whose N basis states is 
used to encode N database items with one-to-one corre¬ 
spondence. Grovers algorithm starts with the state |0) in 
which all qubits are in |0) state. Then it applies Walsh- 
Hadamard transformation W which is nothing but ap¬ 
plication of Hadamard gate on all qubits. After this, we 
get the state which is a uniform superposition of all basis 
states, i.e. 


I 


N-l 


2 = 0 


( 1 ) 


Grovers algorithm then successively iterates the Grover 
operator 


g = Wl^WIr (2) 

on above state to get the basis state \r) encoding the 
solution of search problem. 

Here Iq and A are the selective phase inversions of 
these two states. Mathematically, they are written as 

4 = 12|6)(0| ; A = 12|r)(r| . (3) 

The number of iterations of Q required by Grovers algo¬ 
rithm is (7r/4)-\/]V assuming that there is a unique solu¬ 
tion r. This is quadratically faster than classical search 
algorithms which take 0{N) time steps. Grovers algo¬ 
rithm is proved to be strictly optimal 0. Out of the 
two selective transformations, Jg is easy to implement as 
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we know the state |0). But we dont know the solution 
state |r) in advance and A is implemented using an ora¬ 
cle transformation. Basically, for all |i) in superposition, 
the oracle computes a binary function f{i) which is 1 if 
\i) = |r) else 0. Then the computed value of f{i) is used 
to selectively invert the phase of |i = r). 

In the clause satisfaction problems, we have n Boolean 
variables Xfc (fc G {I, 2,..., n}) and each variable can take 
two values 0 or 1. Let the index i (z G {0,1,..., AI — I}) 
denote the different possible assignments of these vari¬ 
ables where N = 2”. We have m clauses Cj {j G 
{1,2,..., m}) where satisfaction of each clause depends 
upon values of a subset of n variables. Typically this 
subset involves few number of variables. For example, 
in the widely studied NP-complete problem of 3SAT, a 
clause is satisfied if OR of 3 terms is I, where each term 
is either a Boolean variable Xk or its negation Xk (xu = 1 
if = 0 else Xk = 0). The solution is a particular as¬ 
signment i = r which satisfies all m clauses. 

Obviously, in such problems, the function /(z) is AND 
of m different binary functions /j(z) where fj{i) = 1 if 
clause Cj is satisfied else 0. To use Grovers algorithm, we 
need to compute /(z) for all |z) in superposition to im¬ 
plement A- As /(z) is AND of all fj{i)s there has to 
be some coupling between individual computation cir¬ 
cuits corresponding to each These couplings may 

add significantly to physical implementation challenges of 
Grovers algorithm depending upon the kind of hardware 
that will be used for quantum computing in future. 

In this paper, we present an alternative quantum 
search algorithm which is naturally tailored to such kind 
of problems. It relaxes the requirement of coupling com¬ 
putation circuits for each fj{i)s as we dont need to com¬ 
pute AND of all fj{i)s. In next section, we present 
the algorithm and we present its analysis in Section HI. 
The analysis mainly uses the results of general quantum 
search algorithm presented in Q. We then discuss and 
conclude in Section IV. 
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II. ALGORITHM 

Let Ti.N denote the Hilbert space of dimension iV = 2" 
of n qubits where each qubit represents one of the given 
Boolean variables Xk- We attach an ancilla qubit to 
this system and let 'H 2 denote the corresponding two- 
dimensional Hilbert space. We work in the 2N dimen¬ 
sional joint Hilbert space H — 

For each clause j, we define an operator Dj which 
computes the binary function fj{i) corresponding to 
clause and then applies a controlled transformation on 
the ancilla qubit. If fj{i) = 1 then it leaves the ancilla 
qubit unchanged else it applies the single-qubit operator 
Rm on the ancilla qubit, where 

i?m = ea;p(i 7 r/m)| 0 )( 0 | -I- ea;p(—^ 7 r/TO)|l)(l| . (4) 

Obviously, the operator 7)^ acts only on n -|- 1 qubits 
where n is the number of Boolean variables involved in 
clause. The extra one qubit is the ancilla qubit on 
which 'Dj does a controlled operation. The eigenspectrum 
of Dj is of the following form 

2^i(|0)K)) = exp[iTT{l-fj{i))/rn]\0)\i) 

2^i(|l)K)) = ea;p[-i 7 r(l -/j(i))/m]|l)|z) . (5) 

Now consider the operator 

R = Rm'Rm-l ■ ■ ■'R2 'Ri ( 6 ) 

which is basically a product of DjS for all j from 1 to 
TO. As each Dj commute with each other, being diagonal 
in nature, the order of operators doesnt matter for im¬ 
plementing D. Here we note that each unitary operator 
Dj involves computation of fj(i) for only one clause and 
for implementing D, we dont need to couple any of them 
with each other. We just successively apply Dj operators 
for all j in any chosen order. It is easy to see that the 
eigenspectrum of D is given by 

^(|d)K)) = exp[nrui/m]\0)\i) 

T>(|l)|z)) = exp[-nru^/m]\l)\i) . (7) 

Here Ui = ~ is the total number of clauses 

unsatisfied by the assignment i. As Ur = 0 (the solu¬ 
tion satisfies all clauses by definition), we see that D 
has a two-dimensional degenerate eigenspace orthogo¬ 
nally spanned by | 0 )|r) and |l)|r) with eigenvalue 1 . 

Now we present the algorithm. 

(1) Initially, we put all qubits including the ancilla qubit 
in | 0 ) state. 

(2) Then we apply Walsh-Hadamard transform Wn+i on 
all n -I- 1 qubits including the ancilla qubit. We get the 
state |-l-)„+i = |+)|+)n, where |-|-) = l/\/2(|0) -I- |1)) is 
the ancilla qubit state and |-|-)„ = {I/Vn) J2i K) is the 
uniform superposition of all basis states |z) of Hat. 

(3) We perform q times iteration of the operator A = 
Wn+ilQ (jWn+iD on the initial state |-|-)„+i. Here D is 
given by Eq. © and /q g is the selective phase inversion 


of the |0)|0) state. It is easy to see that Wn+il^ ^Wn+i 
is nothing but the selective phase inversion of the 

initial state |-|-)„+i. So we can write the operator A as 

A = I+^n+iR (8) 

After q iterations, we get the state \wq) = A‘^\+)n+i- 

(4) As we show in next section, there exists a q for which 
\wq) is close to the state |-l-)|r) and measuring \wq) will 
output the solution r with a significant probability. 

(5) We repeat the entire algorithm few times to get the 
solution with probability very close to 1. 

III. ANALYSIS 

To prove the working of algorithm, we need to prove 
that A'^\+)n+i is close to |-|-)|r)state for some q. To 
prove this, we prove the reverse operation, i.e. we prove 
that (Ai)'^|-l-)|r) is close to the initial state |-l-)ra+i state. 
From ([8|), we have 

= pt/+,n+l . (9) 

Note that I+,n+i is a self-inverse transformation. The 
above operator is just a special kind of the general quan¬ 
tum search operator S = D^lf that we have analyzed in 
3 . We just need to make the following substitutions in 

5 , 

|s) -)> |-l-)|r), |t) |+)n-|-l, 

5 —>■ (()—>■ TT . (10) 

We refer the readers to original paper Q for the details 
of analysis as here we just use the results of analysis. 

In the analysis of by convention, the initial state 
|s) is the eigenstate of Dg with eigenvalue 1. In case of 
a degenerate eigenspace orthonormally spanned by the 
states 1^), the initial state is chosen according to Eq. 
(1) of 3 3-s 

\s) =-J2{sm\t)\s^) , ^ Ks™|t)p . ( 11 ) 

m m 

Note that here the notation to is used just to conform to 
the analysis of Q. It should not be confused with the 
notation of this paper where to denotes the number of 
clauses. In our case, D has a 2-dimensional degenerate 
eigenspace spanned by | 0 )|r) and 11 ) 1 ^) with eigenvalue 
1. As \t) = |-|-)n-i-i in our case, it is easy to check that 
in (US, the state |s) becomes |-l-)|r) just to satisfy (fTTl) . 
Also, it is easy to check that in our algorithm, 

a = |(r|IT| 6 )| = l/\/]v . (12) 

Using the analysis of Q further, we find that only two 
eigenstates |A±) with the corresponding eigenvalues 
of A^ are relevant for our algorithm as the state |-l-)|r) is 
almost completely spanned by them. These eigenvalues 
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A± are given by Eq. (12) of 0. In our case, this equation 
becomes 


A± = ±— %={tanr])^^ cot 2r] = ■ (13) 

where 

B = y/l + A2 , Ap= ^ |(£|+)„+ipcotP y . (14) 

Here \ tj denote the eigenstates of V with corresponding 
eigenvalues exp{t9i). As the sum is over |£) ^ IV') k) 
(here |V') is any state of ancilla qubit) and 0e=^^r = 0 by 
convention, we need to consider only those eigenstates 
with non-zero eigenphases 9e to evaluate this sum. Using 
Eq. d?]) , we find that such eigenstates with corresponding 
eigenphases are 

|0)|i^r) , exp[nTUi/m] 

|l)|i^r) , exp[—n:Ui/m] . (15) 

So we have \{i\+)n+i\^ = 1/2N for all ^ and we get 


Ap = 


1 




2m 


(16) 


It is easy to check that Ai vanishes because the contribu¬ 
tions from eigenstates | 0 )|i ^ r) identically cancels that 
from eigenstates |l)|i ^ r). Also, if denotes the total 
number of assignments i which dont satisfy u out of m 
clauses, then we can write A 2 as 


A2 


I rri 

-^iV„COt2 


7TU 

2m 


(17) 


With Ai = 0, (fT^ indicates that V = j and so 


A± = ± 


2 

bVn ' 


(18) 


With ?7 = 7 r /4 and (j) = tt, Eq. (23) and (24) of 0 gives 
us the initial state |s) and the effect of iterating S on 
|s) in terms of two relevant eigenstates |A±). With our 
substitutions (??), we get 


|+, r) = -^/V2[e^^+/^\X+) - e*^-/2|A_)], (19) 


and 

(Alt)«|+,r) =-*/v/2[e*«'^+|A+)-e*«'^-|A_)], (20) 

where q' = q + ^. 

For q = q^ Ki 7 r/ 2 |A±| = nBy/N/A, the state 
(ylt) 9 m|-|-^ r) is very close to the state given by 


l/v^(|A+) + |A_)). (21) 


As shown in Q, the above state has an amplitude of 
1/B with the target state \t), which in our case is the 
state |-|-)ra+i- So we have proved that 


K+U+i(Alt)-^'^/4|+,r)| = \{+,r\A-B^/^\+)n+^\ = 1/B. 

( 22 ) 

Thus one running of our algorithm succeeds in finding the 
solution with a probability oil/B^ and O(B^) times run¬ 
ning of our algorithm will give the solution with a proba¬ 
bility approaching 1. So we need a total of 0{'K\fNB^/4) 
iterations of A to solve the search problem. 

In general, for clause satisfaction problems B = 0(1). 

For example, in the case of 3SAT problem, we know that 
with probability very close to 1 , a randomly picked as¬ 
signment i will not satisfy m/ 8 ± 0(i/m) clauses. So the 
sum in (IT71) is approximately cot ~ 4.96 = 0 ( 1 ). 

We point out that our analysis of Q holds only if 
|A±| <C 0min, where 0min is the minimum eigenphase of 
T) different from 0. In our case, 0min = tt/to and as we 
found |A±| = 0(l/vVV), this assumption is satisfied as 
long as iV ^ m? which is true in almost all practical 
situations. 


IV. DISCUSSION AND CONCLUSION 


We have presented a quantum search algorithm which 
naturally relates to the clause satisfaction problems. This 
algorithm allows us to find the solution using Oip/N) or¬ 
acle queries without any necessity of coupling the indi¬ 
vidual computation circuits corresponding to individual 
clause satisfactions. This algorithm has a potential to 
save the computational resources required to implement 
the oracle transformation 1^ in Grover's search algorithm. 

In spirit, our a^orithm is quite similar to that pre¬ 
sented by Kato [ 6 |. There Kato has shown that im¬ 
plementation of /g, the selective phase inversion of the 
| 0 ) (all qubits in | 0 )) state is not necessary for quan¬ 
tum search algorithms and this operator can be replaced 
by an operator made up of only single-qubit gates and 
so physically easier to implement. Here we have shown 
that implementation of the selective phase inversion of 
the solution state A, which is physically harder to imple¬ 
ment, is not necessary for quantum search algorithms for 
clause satisfaction problems. Rather, it can be replaced 
by physically easier to implement operators. 

We believe that similar ideas can be used for other 
kind of search problems also to design quantum search 
algorithms which are physically easier to implement. 
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